Submitted, 2011-09-27; Accepted, 2012-02-14 
Preprint typeset using I^T^X style cmulatcapj v. 5/2/11 
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ABSTRACT 

We present observations of the recently discovered comet-like main-belt object P/2010 R2 (La Sagra) 
obtained by Pan-STARRS 1 and the Faulkes Telescope-North on Haleakala in Hawaii, the University 
of Hawaii 2.2 m, Gemini-North, and Keck I telescopes on Mauna Kea, the Danish 1.54 m telescope 
(operated by the MINDSTEp consortium) at La Silla, and the Isaac Newton Telescope on La Palma. 
An antisolar dust tail is observed to be present from August 2010 through February 2011, while a dust 
trail aligned with the object's orbit plane is also observed from December 2010 through August 2011. 
Assuming typical phase darkening behavior, P/La Sagra is seen to increase in brightness by > 1 mag 
between August 2010 and December 2010, suggesting that dust production is ongoing over this period. 
These results strongly suggest that the observed activity is cometary in nature (i.e., driven by the 
sublimation of volatile material), and that P/La Sagra is therefore the most recent main-belt comet 
to be discovered. We find an approximate absolute magnitude for the nucleus of Hr = 17.9 ±0.2 mag, 
corresponding to a nucleus radius of ~ 0.7 km, assuming an albedo of p = 0.05. Comparing the 
observed scattering surface areas of the dust coma to that of the nucleus when P/La Sagra was active, 
we find dust-to-nucleus area ratios of Ad/ An = 30—60, comparable to those computed for fellow main- 
belt comets 238P/Read and P/2008 Rl (Garradd), and one to two orders of magnitude larger than 
for two other main-belt comets (133P/Elst-Pizarro and 176P/LINEAR). Using optical spectroscopy 
to search for CN emission, we do not detect any conclusive evidence of sublimation products (i.e., gas 
emission), finding an upper limit CN production rate of (Jcn < 6 x 10 23 mol s~ x , from which we infer 
an H2O production rate of Qh 2 o < 10 26 mol s _1 . Numerical simulations indicate that P/La Sagra is 
dynamically stable for > 100 Myr, suggesting that it is likely native to its current location and that 
its composition is likely representative of other objects in the same region of the main belt, though the 
relatively close proximity of the 13:6 mean-motion resonance with Jupiter and the (3,-2,-1) three-body 
mean-motion resonance with Jupiter and Saturn mean that dynamical instability on larger timescales 
cannot be ruled out. 
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1. INTRODUCTION 

Comet P/2010 R2 (La Sagra), hereafter P/La Sagra, 
was discovered by J. Nomen in images obtained on 2010 
September 14.9 (UT) u sing the 0.45 m La S agra Observa- 
tory in southern Spain (Nomen et al. 2010 ). With a semi- 
major axis of a — 3.099 AU, eccentricity of e = 0.154, 
inclination of i = 21.39°, and a Tisserand parameter 
(with respect to Jupiter) of Tj = 3.099, it was immedi- 
ately suspected to be a member o f the recently identifie d 
class of main-belt comets (MBCs; Hsieh fc Jewitt| |2006), 
which are objects that exhibit cometary activity likely 
due to the sublimation of volatile ice, yet are dynam- 
ically indistinguishable from main-belt asteroids. Prior 
to the discovery of P/La Sagra, four MBCs — 133P/Elst- 
Pizarro, 176P/LINEAR, 238P/Read, and P/2008 Rl 
(Garradd) — were recognized, while a sixth possible 
MBC (P/2006 VW 139) has just recently been identified 



( Hsieh et al."|2011c ). Two other comet-like main-belt ob- 
jects— P/2010 A2 (LINEAR) and (596) Scheila — have 
also been observed, but as their comet-like morpholo- 
gies are believed to be due to impact-generated ejecta 
clouds (|Jewitt et al.||2010| |2011| |Snodgrass et al.||2010' 



Bodewits et al.||2011|lYang fc Hsien||2011| |Hainaut et al. 
2012p , and not cometary (i.e., sublimation-driven) dust 
emission, they are better characterized as disrupted as- 
teroids. 

Observationally distinguishing MBCs and disrupted 
asteroids is not straightforward as both types of objects 
exhibit visible dust emission in the form of comet-like 
features like comae and tails. Distinctive morphological 
features helped betra y the true natures of P/2010 A2 

and references 



2012a 



and Scheila, however ( Hsieh et al. 
within). In the case 01 P/2010 A'l, a gap between what 
appeared to be the nucleus of the "comet" and the dust 
tail indicated that the observed dust emission was possi- 
bly an impulsively-generated ejecta cloud that produced 
by an impact onto P/2010 A2's surface and was then 
drifting away. This hypothesis was later corroborated by 
numerical dust modeling and other physic al arguments 



(Jewitt et al 
of Scheila 



2010 



2010). In the case 



Snodgrass et al. 
no gap was seen between the nucleus and 
its emitted dust, but instead of forming a single tail, 
the dust emitted from Schei la appeared to form multiple 
plumes. Hsieh et al. (2012a) considered various scenarios 
that could produce multiple dust plumes, finding scenar- 
ios for sublimation-driven dust emission to be physically 
implausible as they required multiple active sites, which 
were argued to be unlikely on Scheila. Instead, a sim- 
pler scenario in which a single oblique impact caused 
dust emission in the form of a hollow cone which was 
then pushed back in t he antisolar direction by radiation 
pressure was favored ( jlshiguro et al. 2011a|b ). In this 
scenario, what appeared to be multiple plumes of dust 
were proposed to simply be optical projection effects (i.e., 
limb brightening of the hollow dust structure). 

Early observation reports for P/La Sagra indicated 
that it exhibited a largely class ical cometary morph ology, 
with a visible coma and tail dNome n et al.||2010|), sug - 



gesting that it could be a true comet (|Hsieh et al.|^012a ), 
and was therefore likely to be just the fifth MB C t6~be 
discovered a t that time. Numerical modeling by |Moreno 
et al. (20111 furthermore showed that the evolution of 
P/La Sagra s dust tail was consistent with sublimation 



driven dust emission, though they did not explicitly test 
impulsive impact-driven emission scenarios. In this pa- 
per, we seek to further assess whether P/La Sagra's dust 
emission is likely to be cometary (sublimation-driven) in 
nature, or if the object could instead be a disrupted as- 
teroid. We also conduct observational and dynamical 
characterization analyses of the object to place it in the 
context of the small but growing known population of 
comet-like objects in the main asteroid belt. 

2. OBSERVATIONS 

Imaging observations of P /La Sagra were made in pho- 
tometric conditions on multiple occasions in 2010 and 
2011 (Table ^ by the 2.0 m Faulkes Telescope North 
(FTN) on Haleakala in Hawaii, the University of Hawaii 
(UH) 2.2 m telescope, the 8 m Gemini-North Observa- 
tory, and the 10 m Keck I telescope, all on Mauna Kea in 
Hawaii, the 1.54 m Danish telescope at La Silla Obser- 
vatory in Chile, and the 2.5 m Isaac Newton Telescope 
(INT) at the Roque de los Muchachos Observatory on 
La Palma in the Canary Islands (Spain). In addition, 
precovery observations of P /La Sagra were also obtained 
with the 1.8 m Pan-STARRS 1 (PS1) survey telescope on 
Haleakala. To ensure proper calibration of PS1 data, we 
also obtained calibration observations using the Lowell 
Observatory 31" telescope in Flagstaff, Arizona. 

PS1 is a wide- field synoptic survey telescope that em- 
ploys a 1.4 gigapixel camera consisting of a mosaic of 
60 orthogonal transfer arrays (OTAs), each consisting of 
64 590 x 598 pixel CCDs, giving a total field of view 
3.2 degrees in diameter (image scale of 0'.'26 pixel -1 ). 
PS1 observations were made using a Sloan Digi tal Sky 
Survey (SP SS) r' -band-like filter designated rpi ( Stubbs| 



et al. 2010). Gemini observations (obtained as part of 



program GN-2011B-Q-17) were mad e using the Gemini 



2004) 



Hook et al 
hree 2048 x 2048 
band filter, and has 



Multi-Object Spectrograph (CMOS; 
in imaging mode. GMOS employs 
EEV CCDs behind a SDSS-like r 
an image scale of / .'1456 pixel" 1 when using 2x2 bin- 
ning as we did for our observations of P/La Sagra. Keck 
observations were made using the Low Re solution Imag- 
ing Spectrometer (LRIS; Oke et aL| |T995) in both imag- 
ing and spectroscopic mode. LRIS employs a Tektronix 
2048 x 2048 CCD with an image scale of Cf/135 pixel" 1 
and standard Kron-Cousins -B-band and i?-band filters. 
The Danish 1.54 m telescope is currently exclusively used 
by the Microlensing Network for the Detection of Small 
Terrestrial Exoplanets (MiNDSTEp) consortium. Obser- 
vations on the Danish 1.54 m were made using the Dan- 
ish Faint Object Spectrograph and Camera (DFOSC) 
focal-reducing imager behind a standard Kron-Cousins 
i?-band filter. The DFOSC imager has pixel dimensions 
of 2102 x 2102 and a pixel scale of 0'.'39 pixel" 1 . FTN 
observations were made using a 1024 x 1024 e2v CCD 
with an image scale of 0V2785 pixel" 1 behind a Bessell 
i?-band filter. 

Except for observations with PS1, FTN, and the INT, 
all other observations were conducted using non-sidereal 
tracking at the apparent rate and direction of motion of 
P/La Sagra on the sky. Sidereal tracking was used by 
PS1 because this was the standard mode of operation 
for survey operations. FTN observations were obtained 
using its web-based remote control interface, through 
which non-sidereal tracking is not possible. The INT 
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is normally capable of tracking solar system objects non- 
sidereally, but was unable to do so during our observing 
run due to mechanical problems. 

We performed standard image preparation (bias sub- 
traction and flat-field reduction) for UH 2.2 m, Gem- 
ini, Keck, Danish 1.54 m, and INT data using Image 
Reduction and Analysis Facility (IRAF) software. For 
Keck data, flat fields were constructed from images of 
the illuminated interior of the Keck I dome, while for 
the other telescopes, dithered images of the twilight sky 
were used to construct flat fields. PS1 data were re duced 
using the s ystem's Image Processing Pipeline (IPP; Mag- 



mer 



2006) and then calibrated using field stars imaged 



later by the Lowell 31" telescope. Calibration of Gem- 
ini data was performed using field star magn itudes pro- 
vided by the Sloan Dig ital Sky Survey (SDSS;|York et al. 



20001. Photometry of [Lando lt (1992) standard stars and 



field stars was performed tor all data using IRAF and 
obtained by measuring net fluxes within circular aper- 
tures, with background sampled from surrounding cir- 
cular annuli. Conversion of magnitudes measured from 
PS1 and Gemini data obtained using SDSS-like r'-band 
filters to their Kron-Cousins i?-band equivalents was ac- 
complished using the transformation equations derived 
by R. Lupton and made available on the SDSS website 
(http://www.sdss.org/). Comet photometry was per- 
formed using both rectangular and circular apertures as 
described below, where to avoid dust contamination from 
the comet itself, background sky statistics are measured 
manually in regions of blank sky near, but not adjacent, 
to the object. Several (5 to 10) field stars in the comet 
images were also measured to correct for any extinction 
variation (typically negligible) during each night. 

In addition to imaging observations, we also secured 
optical spectra of P/La Sagra on 2010 October 05 with 
LRIS on Keck. We adopted LRIS's l'.'0-wide long-slit 
mask, the 400/3400 grism, and the 460 dichroic on the 
blue side, giving a dispersion of 1.08 A pixel -1 and a 
spectral resolution of approximately 7 A. A total of 3600 s 
of data were obtained from four integrations on the comet 
while it was at an airmass of ~ 1.1. Unfortunately, when 
taking calibration images, our observations were inter- 
rupted due to fog and we were forced to close the dome 
for over three hours. As such, we were only able to record 
a spectrum of a flux standard star and failed to obtain 
any spectra of nearby G-type standard stars, complicat- 
ing our reduction and analysis (Section 3.3.1). Data re- 
duction was accomplished using IRAF. 

3. RESULTS & ANALYSIS 

3.1. Morphological Analysis 

For each night during which observations were ob- 
tained, we construct composite images by aligning each 
night's set of images on the object's photocenter and 
adding them together. We see that P/La Sagra is 
clearly morphologically cometary between August 2010 
and February 2011 (Figure [2]), with a dust tail extending 
up to ~ 30" (~4x 10 4 km at the geocentric distance of 
the object at the time) from the nucleus's photocenter. 
The projected orientation of this dust tail on the plane 
of the sky is seen to rotate counterclockwise throughout 
this period, lagging slightly behind the similarly rotating 
antisolar vector (as projected on the sky). From August 



2010 through November 2010, we observe a single dust 
tail aligned with the projection of the antisolar vector 
in the sky, while between December 2010 and February 
2011, an apparent dust trail aligned with the projection 
of the heliocentric velocity vector (i.e., the plane of the 
object's orbit) appears in addition to the antisolar dust 
tail. By August 2011 (Figure [3]), active dust emission 
appears to have ceased, though a residual dust trail (pre- 
sumably consisting of large dust particles ejected during 
earlier activity) aligned with the projection of the ob- 
ject's heliocentric velocity vector in the sky is observed 
extending at least ~1 arcmin (~ 1.5 x 10 5 km at the 
geocentric distance of the object at the time) from the 
nucleus. This dust trail peaks in brightness well away 

(~ 20") from the nucleus. 

In the formulation o f dust dynamics outlined by |Fin- 



son & Probstein ( 1968 ), the trajectory of an emitted dust 
particle depends on the orbit of the emitting comet, the 
particle's ejection time, its ejection velocity, and the ra- 
tio, (3, of solar radiation pressure, F ra< i, to solar grav- 
ity, F grav , on the particle. In cases where (3 is small 
(i.e., for massive particles), radiation pressure is negli- 
gible, causing particles follow orbits very close that of 
the comet. Over time, these large grains will spread out 
along the comet's orbit to form a potentially observable 
trail. In cases where f3 is large (i.e., for small particles), 
radiation pressure is dominant over solar gravity, giv- 
ing particles nearly directly antisolar motion in the rest 
frame of the comet. In the case of P/La Sagra, we ob- 
serve both limiting cases: a tail consistently oriented in 
the antisolar direction (likely composed of small, high- 
j3 particles) from August 2010 through February 2011 
(Figure [2), and a trail confined to the orbit plane (likely 
composed of large, low-/3 particles) from December 2010 
through August 2011 (Figures [2] and [3]). This assessment 
is in agreement wit h the modeling results obtained by 



Moreno et al. (2011) who found dust grain sizes ranging 
from ~5 /Ltm to ~1 cm. Similar dust structures compris- 
ing widely- varying pa rticle sizes have also been observed 
for comets 2P/Encke (|Reach et al.|2000HLisse et al.|2004 1 
and C/ Austin 1990 V flLisse et al ||1998f ~ 

Due to the ease by which they are accelerated by radia- 
tion pressure, the small particles comprising an antisolar 
tail are expected to dissipate quickly (on timesc ales of 
several wee ks at the distance of the main belt; cf. Hsieh 
et al. 2004| unless they are continually replenished. No 



such separation of the antisolar tail from the nucleus is 
observed at any time from August 2010 through Febru- 
ary 2011, a span of 5.5 months, strongly suggesting that 
continuous dust production was ongoing over this time 
period, consistent with cometary activity. We note that 
the persistence of the orbit-aligned dust trail until Au- 
gust 2011 does not mean that dust production contin- 
ued up until this time as well. In this latter case, the 
longevity of the dust trail is instead likely due to the large 
sizes and therefore slow dis sipation rates of the particles 
in the trail (cf. P/2010 A2; |Jewitt et al.|2010| [Snodgrass 
eTalpOlOl l, not continuous replenishment 



Qualitatively, we find that P /La Sagra's morphological 
appearance and evolution strongly suggest that its dust 
emission is cometary in nature. We caution however that 
more detailed dust modeling (similar that performed by 



Moreno et al.|20TH but also considering impulsive emis- 
sion scenarios) will be required for quantitative confir- 
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mation of the conclusions reached here. 

3.2. Photometric Analysis 

3.2.1. Nucleus Size, Shape, and Color 

Vigorous activity, such as that exhibited by P/La 
Sagra during the majority of our 2010 and 2011 observa- 
tions, normally interferes with the measurement of phys- 
ical properties for cometary nuclei (e.g., size, rotation pe- 
riod, and phase darkening behavior) due to obscuration 
caused by coma when such activity is present. As such, 
in order to accurately measure the physical properties 
of comet nuclei, it is typically necessary to observe them 
once activity has mostly or entirely ceased, as we did with 
P/La Sagra in August 2011 using the Keck I and Gem- 
ini telescopes (Table [2]) . We measure a mean apparent 
i?-band magnitude (averaged in flux space and weighted 
by nominal uncertainties) of m = 23.9 ± 0.1 mag, with 
a magnitude variation between the two nights of Am = 
0.7 ± 0.1 mag (m avg = 24.1 ± 0.1 mag on 2011 August 
26, and m avg = 23.4 ± 0.1 mag on 2011 Aug 31). Pro- 
vided that the photometric variation between these two 
nights is real (intra-night variations are on the order of 
~0.2 mag, with observations on each night only spanning 
30 min or less) and due to rotation, the minimum ratio 
of the axes of the body as projected in the plane of the 
sky is then given by 



0.4Am 



(1) 



giving us a lower limit (due to unknown projection ef- 
fects) to the axis ratio for P/La Sagra of a-N/bN = 
1.9 ± 0.3. Given the extremely limited data set on which 
this calculation is based, this axis ratio should of course 
be considered as a preliminary result only. Better con- 
straints on the shape and structure of the nucleus of P/La 
Sagra will require additional observational sampling of 
the object's lightcurve (likely requiring large telescopes, 
given the small size of the object), ideally at multiple 
observing geometries. 

Using the range of slope parameter (G) values in the 
IAU H , G photometric system computed for fellow MBCs 
133P, 176P, 238P, and P/Garradd (-0 03 < G R < 0.26 



Hsich et al. 2009a 2011b Maclennan & Hsieh 



2012), 



we find an estimated absolute magnitude of 
17.9 ± 0.2 mag. Adopting a i?-ban d albedo of pr = 0.0 5 



(as measured for 133P and 176P; Hsieh et al 
we then find an approximate effective nucleus radius of 
rjv = 0.7±0.1 km, making it comparable in size to 238P, 
about three times larger than P/Garradd, and about 
three times smaller than 133P and 176P (Table]!]). Using 
the approximate axis ratio derived above and assuming a 
triaxial ellipsoidal shape, we find semimajor and semimi- 
nor axes for the nucleus of a at = 1.0 km and — 0.5 km, 
respectively. 

We note that photometric variations on the scale of 
Am = 0.7 mag were not seen at any time while P/La 
Sagra was active (August 2010 through February 2011) 
for any observations that were comparably closely spaced 
in time as the August 2011 observations analyzed above. 
For example, observations on four separate occasions be- 
tween 2010 Sep 08 and Sep 29 exhibited a maximum 
variation of just ~ 0.1 mag (Table [3]). This lack of cor- 
roborating detections of significant photometric variation 
during that period does not invalidate our detection of a 



Am = 0.7 mag variation in August 2011, however, since 
rotational variations for a nucleus obscured by a steady- 
state coma (assumed to be optically thin) will be effec- 
tively damped by that coma. The extent of this damping 
can be determined using 



Am obs = 2.5 log 



F 



obs,max 



F, 



obs 



= 2.5 log 



TV, max 



F. 



N.i 



F H 



(2) 

where Am b s is the observed photometric range, -F bs,max 
and -Fobs, min are t ne maximum and minimum fluxes ob- 
served for the active comet, i*jv,max an d FjV^min are the 
maximum and minimum fluxes for which the nucleus is 
responsible, and the flux due to dust, Fd, is assumed to 
be constant ( Hsieh et al.||2011a ). 

Assuming that the nucleus and dustgrains have sim- 
ilar albedos, flux terms in Equation [2] can be consid- 
ered equivalent to the corresponding scattering surface 
areas of the comet components be ing co nsidered. We 
find in the following section (Section |3.2.2 ) that the ratio 
of scattering surface areas of dust and the nucleus dur- 
ing the four September observations discussed above was 
Ai/A/v ~ 30 (Table [3]). Thus, using A^ max = ira^r^, 
Aw,min = nb N r N , and A d = irr 2 N , we find Am obs = 
0.025 mag. In other words, during this period, even the 
minor ~ 0.1 mag photometric variation we observed for 
the object is unlikely to be due to the rotation of the 
nucleus. Instead, we surmise that the observed variation 
is more likely due to fluctuations in the seeing conditions 
between the different nights, causing fluctuations in the 
amount of coma contained within the fixed 5'.'0 photom- 
etry apertures we used to measure the near-nucleus flux 
of the comet in all of our images. We also conclude that 
our non-detection of significant photometric variations 
in data obtained when P/La Sagra was active is not in- 
consistent with our apparent detection of a significant 
photometric variation while the object was inactive. 

Multi-filter observations using LRIS on Keck I also per- 
mitted us to measure a B — R = 0.9 ± 0.1 mag color for 
P/La Sagra's inactive nucleus, which within our uncer- 
tainties, is consistent with the color of t he Sun (B — R = 



1.03 mag; Hartmann et al. 



1982 1990). The colors we 



measure for P/La Sagra's dust coma while it was active 
are also approximately solar (B — R = 1.02 ± 0.05 mag 
on 2010 October 5 and B - R = 0.94 ± 0.05 mag on 
2009a), 2011 Feb 3), in go od agreement with colors measured for 



the other MB Cs (|Jewitt et al. 
2010b 



2011a). 



2009] |Hsieh et~aT1|2009b[ 



3.2.2. Activity Strength 



Using the nucleus size computed above in Section[3.2.1 



the expected apparent nucleus magnitude for P /La Sagra 
at the time of our first observations in August 2010 is 
mR = 22.2 mag, not = 18.49 mag as measured (Ta- 
ble [3]). This discrepancy indicates that there was signif- 
icant near-nucleus dust contamination even during our 
earliest observations of P/La Sagra. Assuming that the 
phase behavior of P/La Sagra's dust coma is similar to 
that of the nuclei of other MB Cs who se phase functions 
have been measured (cf. Section 3.2.1| since we lack other 
meaningful constraints), we calculate the comet's corre- 
sponding absolute magnitudes, m to t(l, 1,0), at heliocen- 
tric and geocentric distances of R = A = 1 AU and 
phase angles of a — 0° for each set of photometric mea- 
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surements (using photometry apertures with radii of 5'.'0) 
of the nucleus (Table J3]). By comparing these absolute 
magnitudes to our estimated absolute magnitude for the 
nucleus, we then compute the ratio of the scattering sur- 
face area of dust in the coma, Ad, to that of the nucleus, 
An, using 



1 - lo OA { m tot(iA,o)-H R ) 

X00-4(m tot (l,l,0)-H R ) 



(3) 



Then, assuming An ~ 7rr|, and a = 10 /im- r adius grains 
with bulk densities of p = 1300 kg m~ 3 (cf. |Hsieh et ah 
2004 Jutzi et al.|2008 ), we estimate the total dust mass, 
Mdi in the coma using 



(4) 



We list the values we compute for Ad /An and Md using 
the above procedures in Table [3] Comparing the min- 
imum and maximum values of Md during our observa- 
tions, we find that P/La Sagra's coma increases in mass 
at an average net rate of Md ~ 0.1 kg s _1 (likely precise 
to an order of magnitude, at best) between August and 
December 2010. This net increase in total dust mass 
represents the sum of the positive contribution of new 
dust production from the nucleus and the negative con- 
tribution of dust dissipation as individual grains disperse 
beyond the photometry aperture. As such, it represents 
a lower limit to the true dust production rate of the nu- 
cleus. 

In addition to performing nucleus photometry, we also 
measure the total flux from the comet in our compos- 
ite images from each night. We do so by using rectan- 
gular photometry apertures enclosing the entire visible 
dust cloud and oriented to avoid field star contamina- 
tion. Background sky levels were then measured from 
nearby areas of blank sky and subtracted to obtain net 
fluxes, which were then calibrated using standard stars to 
obtain absolute photometry (mtot)- Unfortunately, the 
faintness of the residual dust trail observed on 2011 Au- 
gust 26 and 31, and extensive contamination by nearby 
bright field stars means that we are unable to perform 
these measurements for data obtained on those dates. 
Repeating the mass production rate analysis described 
above, we find that P/La Sagra's coma and tail increase 
in mass at about twice the rate measured for the near- 
nucleus coma, or an average net rate of Md ~ 0.2 kg s _1 . 

Since this analysis does not rely on a fixed aperture 
size, this net mass loss rate is closer to the true dust pro- 
duction rate of the nucleus than that computed above, 
though still represents a lower limit due to the unknown 
dissipation rate of the dust as individual grains disperse 
from the coma and tail and become undetectable above 
the sky background. Comparing with other MBCs, we 
nonetheless find that P/La Sagra has a dust prod uction 
rate compar able to or larger than that of 238P ( |Hsieh| 
et al. 2009b), and at least an order of magnitude large r 



than those of 133P and 176P (Hsieh et al. 112004 2011a 



For reference, we also use our measurements ot the 
near-nucleus flux of the comet (using 5'.'0 photometry 
apertures) in each set of observations to e stimate the dust 
contribution as parameterized by Afp ( |A'Hearn et al. 
1984 1. These results are also shown in Table [3] and again 



indicate significantly stronger activity for P/La Sagra as 
compared to 133P, for which A fp was also systematically 
tabulated ( |Hsieh et aL|2010b[ ). 

We also attempt to reproduce our observed photom- 
etry using a sublim ation model similar to that used by 
( |Hsieh et al.|2011b[ ) to analyze the activity of 238P. This 
effort was hampered, however, by our lack of observa- 
tional constraints on the start time of activity and the 
rate of its eventual decline. We do find that steadily in- 
creasing activity during the post-perihelion portion of 
P/La Sagra's orbit cannot be modeled by surface ice 
sublimation, indicating that the ice must instead reside 
at some depth below the surf ace of the nucleus, s imilar 



to what was found for 238P (Hsieh et al. 2011b). Un 



fortunately, the numerous free parameters that can be 
modified in our model mean that this finding is the only 
meaningful result that can be derived at this time. Bet- 
ter sampling of the activity over a longer time baseline 
will be required to gain additional insight into the na- 
ture of P/La Sagra's activity via this method, and as 
such, should be considered an observational priority dur- 
ing its next expected period of activity (likely beginning 
in mid-2015). 

We note that if P/La Sagra exhibits activity over a 
period of ~ 1 ye ar during each of its perihelion passages 
(similar to 133P; Hsieh et al.|2004 ), assuming an average 
mass loss rate of Md = 0.2 kg s _1 , a total mass loss 
of ~ 10 7 kg per orbit would be expect ed. A ssuming a 
nucleus size of vn = 0.7 km (Section 3.2.1) and p = 
1300 kg m~ 3 , such a mass loss rate should consume the 
entire nucleus over the next ~ 10 5 orbits, or < 10 6 yr. 
However, activity in MBCs (even over just a portion of 



their orbits) is expected to be transient (cf. |Hsieh et al 
2004[ ), p erhaps amou nting to just 1% of an object's total 
lifetime ( Hsieh|2009 ). As such, we expect a timescale for 
nucleus disintegration due to cometary activity closer to 
~ 10 8 yr (Hsieh 2009), or approximately comparable to 
the ~ 10 8 -10 y yr collisional lifetime expected f or a main- 
be lt aste roid the size of P/La Sagra's nucleus (Bottke et 
aTpK)5|). 



3.2.3. Comparison to Other MBCs 

To further understand P/La Sagra's place among the 
population of known MBCs, we plot Ad/AN for the 
first five known MBCs as functions of true anomaly, 
v (Figure [5k). We note a few interesting features of 
this plot. First, we find quantitative confirmation that 
for their nucleus sizes, 238P and P/Garradd exhibit far 
more dust when active than 133P and 176P (also see 
Table [l}. This fact has of course always been qualita- 
tively apparent. The observed morphologies of 238P and 
P /Garr add, when active, resemble those of highly activ e 



comets (Hsieh et al 
whereas the nuclei o 
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Jewitt et al. 2010) 



dominate the sur- 
face brightness profiles of those objects, even at the peaks 
of their active periods. Only comparatively faint dust 
tails and small amounts of photometrically-detected un- 
resolved near-nucleus dust hav e attested to the cometary 



nature s of these latter MBCs ( |Hsieh et al.||2004| |2010b 
2011a). As for P/La Sagra, we find that it exceeds even 
238P and P/Garradd in the amount of dust it produces 
for the size of its nucleus, producing ~ 3 times as much 
dust relative to its nucleus size at its peak as 238P and 



G 



P/Garradd produced at their peaks, which in turn were 
~ 50 times as productive relative to their nucleus sizes 
as 133P and 176P. 

We caution that given the different geocentric dis- 
tances at which the observations in question were made, 
the near-nucleus photometry performed for each MBC 
measures dust contained within different physical dis- 
tances from the nucleus. As such, considering that other 
physical characteristics of the dust emission (e.g., grain 
sizes, bulk grain densities, and ejection velocities) for 
each MBC are unknown and likely also vary from ob- 
ject to object, these comparisons are best interpreted in 
terms of orders of magnitude. We therefore find 238P, 
P/Garradd, and P/La Sagra to have similar levels of 
activity which are 1-2 orders of magnitude higher than 
those of 133P and 176P, implying similar ratios of frac- 
tional active surface areas for these objects. 

To quantify the total amount of ejected dust exhib- 
ited by each comet, we also multiply the Ad /An values 
for each comet by estimates of Am for each object (de- 
rived from nucleus radii listed in Table [T]) to obtain ab- 
solute values for Ad- We plot the results as functions of 
true anomaly in Figure ^5jp. Interestingly, we find that 
despite the high Ad/A^ values measured for 238P and 
P/Garradd, the total absolute amount of dust in those 
comets is approximately comparable to the amounts of 
dust present for 133P and 176P, which visually appear 
far less active in comparison. P/La Sagra is seen to pro- 
duce roughly an order of magnitude more dust in abso- 
lute terms than any of the other MBCs (assuming similar 
particle size distributions). 

We therefore find the interesting result that within an 
order of magnitude, all five MBCs considered here ex- 
hibit similar total amounts of visible dust in their comae 
and tails, despite the wide variation in their observed 
morphologies. Determining whether this finding is phys- 
ically significant will require similar dust mass measure- 
ments for newly-discovered MBCs to investigate whether 
this trend is maintained as the sample of known objects 
grows. We suggest that future research investigating the 
validity and implications of this possible trend may pro- 
vide valuable insights into the origin and evolution of 
MBC activity. 

3.3. Spectroscopic Analysis 
3.3.1. Data Reduction 

In the optical region of the electromagnetic spectrum, 
the emission band of the CN radical at 3880 A is among 
the strongest among cometary gaseous species and there- 
fore is considered to be the most sensitive probe for de- 
tecting gas resulting from sublimating ice. We search 
for evidence of the presence of the CN band in a one 
dimensional spectrum extracted from the spectral image 
using a l'.'O x 7'.'5 rectangular aperture centered on the 
continuum. We use a relatively wide extraction aper- 
ture to incorporate as much light as possible from the 
surrounding coma. The sky background is reconstructed 
and subtracted using flanking regions extending from 9" 
to 15" from the nucleus. 

In Figure [6j we present the spectrum of P/La Sagra in 
the relevant spectral region from - 3700 A to - 4000 A. 
The continuum of the comet consists of scattered sun- 
light from dust particles in the inner coma and shows 



several absorption features in the Solar spectrum, such 
as the prominent Ca H and K absorption lines at 3933 A 
and 3966 A (Figure [6| . The shaded spectral region from 
3830 A to 3905 A indicates where we expect to d etect 
CN emission if any is present ( Cochran et al.|[l992 ) . 

In order to correct for solar spectral features, the spec- 
trum of a G2V solar analog star is needed. As discussed 
in Section [2j we were not able to observe a G2V star 
at the time of our P/La Sagra observations due to poor 
weather conditions. We instead adopt a flux-calibrated 
spectrum of G2V star HD 91163 which was obtained us- 
ing the same instrument but was r ecorded a few month s 
later on UT 2010 December 17 ( |Hsieh et aL||2012a ). 
This spectrum was observed using the 600/4000 grism 
and thus has a higher spectral resolution than our P /La 
Sagra data, however, so the additional step of binning 
the HD 91163 spectrum to the same spectral resolution 
as our P/La Sagra spectrum is required. The scaled spec- 
trum of HD 91163 is shown as the blue dashed line in 
Figure [6] We note that the spectra of P/La Sagra and 
HD 91163 are quite similar, particularly in the shaded 
region. However, minor discrepancies are observed at 
wavelengths A < 3800 A and A > 3980 A. These discrep- 
ancies may be due to the time and airmass differences 
between the comet observations and standard star ob- 
servations. We find no evidence of CN emission in the 
shaded wavelength region. 

3.3.2. Upper Limits on Gas Production Rates 

The continuum-removed spectrum of P/La Sagra is 
shown in Figure [7J In principle, this spectrum shows 
emission from atoms and molecules in the inner coma. 
We measure the standard errors in the residual spectrum 
in three wavelength regions 70 A in width: 3760 A to 
3830 A, 3830 A to 3900 A, and 3900 A to 3970 A. The 
standard errors on either side of the expected CN band 

are 5.1 x 10~ 18 and 2.8 x 10~ 18 erg cm" 2 s" 1 A~ , re- 
spectively, shown as red error bars in Figure [7j The stan- 
dard error in the region of the expected CN band is 4.2 x 

10 -18 erg cm -2 s -1 A . We choose the largest standard 



error of 5.1 x 10 18 erg cm 



A" 1 



as a conservative 



estimate of our observational uncertainty in the region of 
the expected CN emission band. We therefore estimate 

a 3-cr upper limit of 1.5 x 10~ 17 erg cm~ 2 s" 1 A 

We calculate the integrated CN band flux by summing 
the estimated emission flux in the shaded GN region, 
obtaining /cn = 3.7 x 10~ 16 erg cm~ 2 s . We then 
convert the integrated flux to the total number of CN 
molecules using 



4^A 2 / CN 



(5) 
(0 



where A^cn is the number of CN molecules, A and 
R are the geocentric and heliocentric distances, re- 
spectively, and g(R) is the resonance fluorescence ef- 
ficiency, which describes the number of photons scat- 
tered per second per radical for an optically thin coma, 
in erg s _1 molecule -1 at a heliocentric distance of R. 
At the time of our observations, P/La Sagra had a ra- 
dial velocity of R = 1.1 km s _1 , for which g(R — 
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1 AU) 



= 2.86 x 10 
3.7 x KT 16 



13 



(jSchleicher 20101 
r 



Substituting 

f CN = 3.7 x 10- iB erg cm" 2 s" 1 , R = 2.66 AU, and 
A = 1.82 AU, we find an upper limit of iV = 8.3 x 10 25 
CN molecules in the extra ction apertu re. 

A simple Haser model (Haser 19571 was used to de- 
rive the CN production rate based on the number of 
C N molecules. The d etails of the model are described 
in Hsieh et al. (2012a). We find an upper limit to the 
CN production rate of Qcn < 6 x 10 23 mol s" 1 . Taking 
average ratios measured for previously observe d comets 
(logf(WC ?OHl = -2.5; Qoh/Qh ? = 90%) ( [A'Hearn 
et al.|1995 1 , we estimate an upper limit on the water pro- 
duction rate of Qh 2 o < 10 26 mol s^ 1 . Incorporating the 
Afp value calculated fro m dat a taken at the same time 
as these spectra (Section 3.2.2| Table J3|, we find a lower- 



limit dust-to-gas ratio of log[Afp/Qc^} > —22, simi- 
lar to values measure d for other Jupiter-family comets 
(A'Hearn et al. 1995). As this result is only an lower 
limit, however, it is not possible to ascertain from these 
data the degree to which P/La Sagra is depleted in 
volatiles relative to other classical comets from its ap- 
parently long residence time in the inner solar system 
(Section Ktn. 



3.3.3. Estimation of Effective Active Area 

We can estimate the effective active area of the nucleus 
of P /La Sagra by comparing the observed mass loss rate 
with theoretical rates obtained from applying a simple 
thermal model. To calculate the sublimation rate of sur- 
face ice, we solve the energy balance equation for the 
equilibrium temperature on the object's surface: 



S e (l-p) 



J{9) = eoT* + Lyj(T e ) 



(7) 



where 5 Q = 1360 W s _1 is the solar constant, p = 0.05 
is the albedo, R = 2.66 AU, J(9) is a function of the so- 
lar zenith angle 9 which describes the ratio between the 
absorbing and radiating area, e = 0.9 is the emissivity, a 
is the Stefan-Boltzmann constant, L = 2.68 x 10 6 J kg -1 
is the latent heat and T e is the effective surface temper 
ature. The specific sublimation rate, ip (in kg m~ 2 s 
of water is given by 



V> = P(T e 



2nkT P 



(8) 



where m is the mass of a water molecule, k is Boltz- 
man n's constant, and P{T e ) is the saturated vapor pres- 
sure dFanale & Salvail||1984|). 



In the fast rotator approximation (which yields the 
coldest surface temperature), the nucleus is considered 
isothermal and J = (cos 9) — 1/4. The hottest surface 
temperatures are derived in the sub-solar point (i.e., slow 
rotator) approximation, where the absorbing and radiat- 
ing areas are identical and therefore J = 1. We ob- 
tain sublimation rates of 1.3 x 10~ 6 kg m~ 2 s _1 and 
6.1 x 10~ 5 kg m~ 2 s _1 for the cold and hot scenar- 
ios, respectively. Assuming an av erage mass loss rate 
of Md ~ 0.2 kg s _1 (Section 3.2.2), we estimate active 
areas of 1.5 x 10 5 km 2 and 3.3 x 10 3 km 2 , correspond- 
ing to total active surface fractions of ~3% and ^0.05% 
(likely precise to an order of magnitude, at best), for the 
cold and hot cases, assuming dust-to-gas mass ratios of 



unity. These derived active surface fractions are compa- 
ra ble to those of o t her J upiter-family comets measured 
by jA'Hearn et aL (19951. For comparison, both 133P 
and 176 P are estimated to have activ e areas of ~ 10 2 - 
10 3 m 2 ( |Hsieh et al.|2004[ |Hsieh|2009l) corresponding to 
active surface tractions of 0.0002 - 0.002%, several orders 
of magnitude smaller than for P /La Sagra. 

3.4. Dynamical Analysis 
3.4.1. Stability Analysis 

To better understand the likely origin of P/La Sagra, 
we als o consider its dynamical properties. |Jewitt et al.| 
(|2009| raised the possibility that MBC P/Garradd could 



have originated in a different part of the asteroid belt 
from where we see it today. As such, an important issue 
to consider is whether P/La Sagra is dynamically stable 
at its present location, or whether it is possible that it 
may have also originated elsewhere. 

To assess P/La Sagra's dynamical stability, we ran- 
domly generate three sets of 100 test particles, each set 
with Gaussian distributions in orbital element space, cen- 
tered on P/La Sagra's JPL-tabulated osculating orbital 
elements, where each set is characterized by a a value 
equal to the JPL-tabulated uncertainties (as of 2011 Jan- 
uary 1) for each orbital element. To investigate the sta- 
bility of the region (in orbital element space) surrounding 
P /La Sagra, we also generate three additional sets of 100 
Gaussian-distributed test particles each, where each set 
is again centered on P/La Sagra's osculating orbital el- 
ements but is characterized by a a value equal to 100 x 
the JPL-tabulated uncertainties. We then use the N- 
body integration package, Mercury ( Chambers|1999 ), to 
integrate the orbits of each set of test particles forward 
in time for 100 Myr (limitations on available comput- 
ing resources unfortunately prevent us from conducting 
these simulations over significantly longer time periods, 
e.g., 1-2 Gyr, in a reasonable amount of time), where we 
treat the eight major planets as massive particles and all 
test particles are considered to be massless. 

Performing these integrations for all six sets of test 
particles, we find that 7% of the objects in the three 1- 
o test particle sets are ejected (i.e., reach a heliocentric 
distance of > 50 AU) from the asteroid belt over the 
course of 100 Myrs (Figure [8]), and 11% of the objects in 
the three 100-cr test particle sets are ejected (Figure [9]). 
P/La Sagra itself is found to be stable over the 100 Myr 
test period, consistent with its Lyapunov time (Ti yap ), 
estimated by us to be ^200 kyr (comput ed using the pro- 
cedure described in Tsiganis et al.|2003 ) , where a body is 
generally considered to be stable for 'Ji ya p > 100 kyr. The 
l-<7 test particles that are ejected are not localized to par- 
ticular regions of orbital element space near P /La Sagra, 
suggesting that their ejection may be a consequence of 
stochastic encounters with a weak mean-motion reso- 
nance (or multiple resonances), rather than the effect 
of a single dominant resonance. Nearby resonances in- 
clude the "moderate-order" 13:6 mean-motion resonance 
with Jupiter at a = 3.1063 AU and the (3,-2,-1) three- 
body mean-motion resonance (with Jupiter and S aturn) 
at a = 3.0790 AU ( |Nesvorny fc Morbidellil[l998| . Both 
resonances are more than 1-a away in a, given the uncer- 
tainties for P/La Sagra's osculating elements, but could 
conceivably affect the object as its osculating elements 
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fluctuate over time (Figure 10 1. 



In our 100-ct test particle simulations, the effect of 
the (3,-2,-1) three-body resonance is clear, with ~70% 
of the ejected test particles in those simulations having 
initial semimajor axes within 0.01 AU of the resonance. 
The effect of the 13:6 resonance is comparatively less 
pronounced, with only ~20% of the ejected particles in 
these simulations starting within 0.01 AU of the reso- 
nance, despite the fact that this resonance is closer to 
P/La Sagra than the (3,-2,-1) resonance and therefore 
lies in a more highly populated region of our initial test 
particle distribution. P/La Sagra itself has an initial os- 
culating semimajor axis 0.06 AU of the 13:6 resonance, 
and even occasionally crosses the resonance during our 
simulations (Figure 10), yet is found to be stable over 



the 100 Myr test period considered here, suggesting that 
this resonance only weakly affects objects in this region 
of orbital clement space. 

We conclude from this analysis that the region in or- 
bital element space occupied by P/La Sagra is largely 
stable, implying that the object is likely to be native 
to its present-day location. As such, its physical prop- 
erties may be reflective of other objects in this region 
of the asteroid belt. However, between the nearby 13:6 
mean-motion resonance with Jupiter and (3,-2,-1) three- 
body mean-motion resonance with Jupiter and Saturn, 
enough dynamical instability is imparted to the region 
that a non-negligible fraction of simulated test particles 
escape in less than 100 Myr. Due to the non-negligible 
instability of this region, we therefore cannot definitively 
conclude that P/La Sagra formed in situ (and did not 
dynamically evolve onto its present-day orbit from else- 
where in the main belt, or even elsewhere in the so- 
lar system), particularly over timescales longer than the 
100 Myr test period considered in this analysis. We 
further note that our analysis does not consider non- 
gravitational effects such as the Yarkovsky effect or the 
effects of asymmetric mass loss such as that associated 
with cometary jets. As such, in future investigations of 
this object, particularly those which pertain to the com- 
parison of its physical properties to those of other MBCs, 
the possibility that P/La Sagra may not be native to its 
current location cannot be discounted. 

3.4.2. Family Search 

Having investigated the dynamical stability of P/La 
Sagra and its surroundings, we seek to further character- 
ize the object's dynamical environment by searching for 
evidence of an associated asteroid family, i.e., a grouping 
of dynamically linked asteroids likely to be the result of 
either th e catastrophic f ragmentation of a larger parent 
asteroid (| Hirayama|1918 [) , or a significan t cratering even t 
on a still-existmg asteroid (e.g., Vesta; Asphaug 1997). 
Such a search is particularly crucial for this object be- 
cause of studies suggesting that MBCs may preferentially 
be fou nd among families ( Hsieh et al. 2004[ Haghighipour 
2009[ ), a nd possibly specifically among extremely young 
families ( [Nesvorny et al.||2008l |Hsieh||2009j ). The lat- 
ter hypothesis is based on the facts that MBC 133P ap- 
pears to be a member of the < 10-Myr-old Beagle family 
( Nesvorny et al.|[2008 ), and collisional devolatilization of 
near-surface ice of kni-scale MBCs (like 133P and 176P) 
is expected to occur at rates incompatible with the ex- 
istence of present-day activity, unless those MBCs are 



recently-produced fragments o f larger aste roids with sig- 
nificant interior ice reservoirs ( |Hsieh||2009| ). 

Due to P/La Sagra's high inclination, we cannot 
search for an associated dynamical asteroid famil y using 
analytically-determ ined proper orbital elements (Milani 
|fc Knezevic1[T994| , as can be done for objects at lower 
inclinations. The complication stems from the fact that 
analytically-computed proper elements are significantly 
less accurate for objects with high eccentricities or in- 
clinations. The development of numerical methods for 
computing so-called synthetic proper elements for char- 
acterizing the long-term behavior of high-inclin ation and 
high-eccentricit y objects resolved this problem ( Knez evic 
fc Milanj]|2000| ) , an advance that was empl oyed by Gil 
Hutton| ([2006 ) and Novakovic et al. (2011) to conduct 
systematic searches tor families among high-inclination 
asteroids. 

For our analysis of P/La Sagra, we first compute the 
object's own synthetic proper elements — proper semi- 
major axis (cip), proper eccentricity (e p ), and proper 
inclination - - usin g the procedure described in 
Knezevic & Milani|(|2000|) and osculating orbital elements 



from the J -PL Small-Body Database Browser. We list 
these results in Table |4j We then plot these newly com- 
puted synthetic proper elements in a p — e p and a p — i p 
space along with the synthetic elements of all of the fam 
ilies (as well as smaller or less-significant gr oupings des 
ignate d as clumps or clusters) id entified by |Gil-Hutton 
( |2006| ) and |Novakoyic et al.J ( 201 1 ) to check for plausible 
close associations of F/La Sagra w ith any of these previ- 
ously identified groupings (Figure 11). We find no such 
associations. 

Finding no link between P/La Sagra and any currently 
known asteroid families, we turn to our own search for 
dynamically associated asteroids to determine whether a 
clustering specific to P /La Sagra may have missed simply 
because this newly discovered object was not considered 
in previous family searches. To do so, we add our com- 
puted proper elements for P /La Sagra to an updated list 
of sy nthetic proper element s for high-inclination aster 
oids (|Novakovic et al.||2011[), and t hen employ the Hier 



archi cal Clustering Method (HCM; Zappala et al. 1990 



1994 1 to search for dynamically associated asteroids. Re 



suits from this analysis are plotted in Figure [12] where 
we show the number of asteroids associated with P/La 
Sagra as a function of cutoff "distance" (which, despite 
its name, has units of velocity) . 

We find just six objects in P/La Sagra's dynamical 
vicinity (including P/La Sagra itself) within a cutoff dis- 
tance of 120 m s _1 (Table pi). Confirming our search 
(described above) for associations of P/La Sagra with 
previously identified groupings, none of these objects be- 
long to any of the fa milies, clumps, or cluste rs found by 



Gil-Hutton| ()2006| ) or |Novakovic et al. ( |201 1 ) . According 
to the criteria used by Novakovic et al. for the outer belt 
region, this group of objects is neither compact nor large 
enough to be classified as a family or cl ump. At best, it 
can be classified as a cluster, defined by [Novakovi c et al.| 
(20111 as a compact group of objects that could share 



a common physical origin, but does not formally satisfy 
criteria for being a family or clump, usually because the 
number of members that it contains is too small. 
To obtain a quantitative assessment of the likelihood 
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that members of the P/La Sagra group share a common 
physical origin, we use a method employed by Nesvorny] 
& Vokrouhlicky ( 2006 ) to estimate the probability that a 
certain cluster of objects is the result of random chance. 
Confining our analysis to the high-inclination region of 
the outer main belt (2.825 < a p < 3.278 AU, e p < 0.4, 
0.3 < sin(ip) < 0.55), we generate 10,000 random distri- 
butions of 9569 synthetic objects (the same number of ob- 
jects currently known to populate this region) in proper 
orbital element space. We then perform HCM analy- 
sis on these synthetic populations to search for possible 
groupings with at least 6 objects within cutoff distances 
of 95 m s _1 . The results of this analysis show that there 
is a significant probability (~ 30%) that the P/La Sagra 
cluster could in fact be due to chance. While these re- 
sults still leave a ~ 70% chance that this grouping is not 
random and therefore could have a physical origin, more 
information is certainly needed to establish whether or 
not this clustering can be considered real. 

We compute Lyapunov times for members of the P/La 
Sagra group (Table[4]), finding the orbits of all six bodies 
to be stable (7i yap >100 kyr). This fact provides an 
opport unity to apply the bac kward integration method 



2002) to this group. BIM is a 



(BIM; |Nesvorny et al 
means for estimating the age of a <20-Myr-old family 
that is based on the premise that immediately after the 
disruption of a parent body, the orbits of the resulting 
fragments must have been nearly the same, and has been 
successfully applied to est imate the ages of th e Karin 
cluster and Veritas family ( Nesvorny et al.|[2003 ). 

In the case of the group of asteroids surrounding P/La 
Sagra, the purpose of BIM analysis is twofold. If a sus- 
picious clustering in the longitude of the ascending node 
(f2) or the longitude of perihelion (m) can be found, this 
analysis would not only allow us to estimate the age of 
the group, but would also provide support for the hy- 
pothesis that these objects share a physical origin in 
the first place. In the outer belt, however, application 
of this method is generally complicated by large gra- 
dients of secular frequencies with a pi particularly when 
consider ing w in the vicinity of the 2:1 resonance with 
Jupiter ( Nesvorny et al.|2003 ). Although the P/La Sagra 
group is not particularly close to this resonance, the 
gradient dg/da p , where g is the frequency of the lon- 
gitude of perihelion, is still relatively large in this region 
(~0.1° yr _1 AU -1 ). As such, for the purpose of our 
study we only consider f2. 

To investigate whether or not a clustering in fl ex- 
ists, we numerically integrate the orbits of all six known 
members of the P/La Sagra group for 20 Myr backward 
in time. We then calculate the average difference, (Af2) , 
calculated for all possible combinations of objects, where 
AO = tti — Clj, and i and j denote two bodies under 
consideration. Plotting the results of this analysis (Fig- 
ure fl3 ) , we find three clusterings within 60 degrees. The 



clusterings at about 4.9 and 18.2 Myr ago are less sig- 
nificant than the one occurring about 6.2 Myr ago. This 
implies that the age of the cluster is likely close to 6.2 
Myr, but due to the small number of asteroids involved, 
it is not possible to draw a definitive conclusion about 
the physical significance of this cluster. 

The results shown in Figure [13] are obtained in a purely 
gravitational model. Including other forces (e.g., the 
Yarkovsky effect) could affect the positions and depths of 



various minima (cf. Novakovic 2010[ ). We believe, how- 
ever, that such an analysis should be deferred pending 
the discovery of more members of the cluster so that more 
statistically meaningful results can be obtained than are 
currently achievable with just six possible cluster mem- 
bers. In the meantime, the determination of colors or 
acquisition of spectra of the currently suspected cluster 
members could provide some preliminary insights into 
the likelihood that these objects are physically related. 

4. DISCUSSION 

To date, the only published an alysis of P /La Sagr a 
prior to this work was presented by Moreno et al. (2011 ), 
who focused on numerical modeling of object's dust emis- 
sion. They found that the activity can be plausibly mod- 
eled by an episode of continuous dust emission, consis- 
tent with a sublimation-driven dust ejection process and 
inconsistent with impulsive dust ejection as would be ex- 
pected from an impact, beginning shortly before the ob- 
ject's perihelion passage and persisting over >7 months. 
Furthermore, they found that the dust cloud's morphol- 
ogy was best reproduced by anistropic dust emission 
and indicated that their modeling suggests that P/La 
Sagra's activity could be seasonally mo dulated in a man- 
ner sim ilar to that of other MBCs (e.g., |Hsieh et al.|2004 



2011a). In ou r photometric analysis of the dust clouc 
(Section 3.2.2), we found steadily increasing dust mass 
between August 2011 and December 2011, indicative of 
ongoing d ust production , and consistent with the con- 
clusion by Moreno et al. (2011) that P/La Sagra's dust 
emission appears to have persisted over an extended pe- 
riod of time, and therefore was likely to be cometary (i.e., 
sublimation-driven) in nature. 

While the conclusion that P/La Sagra's dust emission 
is driven by the sublimation of volatile ices appears to 
be at od ds w ith the results of our spectroscopic analysis 
(Section 3.3), we note that the extreme difficulty of de- 
tecting sublimation products in a cometary coma at the 
distance of the main belt means that the non-detection 
of gas emission in such observations cannot be consid- 
ered conclusive evi dence against the presence of gas (cf. 
Hsieh et al.||2012a |. Rather, such non-detections should 
only be interpreted as evidence that the levels of gas 
emission for all MBCs for which spectroscopy has b een 
obtained to date (133P, P/Garradd, and P/La Sagra; Li- 
candro et al.|2011 Jewitt et al.j2009 and this work) have 
simply been below the sensitivity levels of those observa- 
tions. However, the aforementioned observations were all 
obtained using large ground-based telescopes (the Euro- 
pean Southern Observatory's 8 m Very Large Telescope 
in Chile, and the 10 m Keck I telescope in Hawaii), and so 
future attempts to detect gas in MBCs may need to focus 
on significantly more active MBCs (i.e., with higher gas 
emission levels) since significantly more powerful observ- 
ing facilities are not currently available. It may also be 
possible that given the large distances and low activity 
levels (relative to other comets observed at much smaller 
heliocentric and geocentric distances) of the MBCs, con- 
clusive evidence of gas emission may only ever be ob- 
tainable from a close-approaching spacecraft equipped 
to make such measurements. 
D espit e the lack of a successful gas detection (Sec- 



tion 



3.3), we conclud e from the numerical modeling of 
Moreno et al.| (2011) and the morphological and pho- 
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tom ctric analyses presented in this work (Sections 3.1 
and 3.2) that P/La Sagra's activity is likely to be 
sublimation-driven. As such, the object is likely to be 
a bona fide MBC like 133P, 176P, 238P, and P/Garradd, 
where dust is ejected via gas drag from sublimating ice, 
and not as ejecta from a recent impact as is suspected 
in the cases o f disrupted asteroids, P/2010 A2 and (596) 



Scheila (e.g., Jewitt et al. 



2010 



Bodewits et al.||2011). 



2010 2011 Snodgrass et al. 



Definitive criteria for identifying the true nature of 
comet-like activity remain elusive, howev er (cf. Hsieh et 



al.|| 2012a[ ). In the case of P/La Sagra, Moreno et al 
] 2T)lip show that comet-like dust emission can plausi- 
bly explain the appearance of P/La Sagra's dust tail, 
but by no means show that their solution is unique, a 
common weakness in virtually all dust modeling efforts 
where models are underconstrained by available observa- 
tions. Given the case of P/2010 A2's long -lived dust tail 
( |Jewitt et al.||2010[ [Snodgrass et aL|2010[ ), the longevity 
of P/La Sagra's activity is likewise not incontrovertible 
evidence of sublimation-driven dust emission. Even the 
apparently steadily inc reasin g dust mass in P /La Sagra's 
coma and tail (Section |3. 2. 2 1 could be considered incon- 
clusive evidence of sublimation given the potential that 
fragmenting dust particles ejected in an impulsive event 
could mimic ongoing dust production as in the case of th e 
2007 outburst of comet 17P/Holmes ( |Hsieh et al |2010a[ ). 

At this time, recurrent activity, separated by interven- 
ing periods of inactivity, appears to be the only reliably 
observable indicator of sublimation-driven dust emission 
that cannot also be plausibly explai ned by a scenario in - 



volving impact-driven dust ejection ( Hsieh et al. [2011c 
Repeated impacts on individual asteroids are unexpected 
from either theoretical calculations or empirical obser- 
vations. The tendency of episodes of repeated activity 
to occur over similar arcs of each object's orbit makes 
such behavior even more difficult to explain as the ac- 
tion of random impacts rather than the effect of periodic 
cometary sublimation. Two of the six currently known 
MBCs have been observed to exhibit recurrent activ- 
ity to date (133P and 238P; 



2011b), while 176P is expected to exhibit renewed activ 



Hsieh et al. 2004 



2010b 



ity shortly ( Hsieh et al.|2011a ). Confirmation of renewed 



activity in 176F will be crucial for supporting the hypoth- 
esis that its activity is cometary in nature. The remain- 
ing three MBCs (P/Garradd, P/La Sagra, and P/2006 
VWj.39), on the other hand, were actually discovered so 
recently that they have not yet completed full orbits since 
being discovered, but as they do (P/Garradd in early 
2013, P/La Sagra in mid-2015, and P/2006 VW139 in 
mid-2016), observations to search for recurrent activity 
are highly encouraged. 

5. SUMMARY & CONCLUSIONS 

We present a wide-ranging study of the recently discov- 
ered comet P/2010 R2 (La Sagra) including photometric 
and spectroscopic observations and dynamical analyses 
of the object. Key results are as follows: 

1. A year-long observational monitoring campaign of 
P/La Sagra from August 2010 to August 2011 show 
that the morphology of its dust cloud undergoes 
substantial evolution, exhibiting a strong coma and 
a single tail from August 2010 to November 2010, 



and a coma and two tails from December 2010 to 
February 2011, before diminishing to a faint, de- 
tached dust trail by August 2011. A long-lived an- 
tisolar dust tail and the steadily increasing scatter- 
ing surface area measured for the dust cloud over 
much of this period suggests that dust production 
was ongoing over a period of several months, be- 
havior consistent with dust ejection via the subli- 
mation of volatile ice and inconsistent with dust 
ejection via the impulsive action of an impact on 
the surface of P /La Sagra. 

2. We estimate that the nucleus of P/La Sagra has 
an absolute magnitude of H R — 17.9 ± 0.3 mag, 
corresponding to an approximate effective radius 
of r N = 0.7 ± 0.1 km. The B - R color of the 
nucleus {B — R = 0.9 ± 0.1 mag) is approximately 
solar, as are the B — R colors of the active comet 
measured on two separate occasions in 2010. 

3. Optical spectroscopic observations of the active 
P/La Sagra 20 days after it was discovered reveal 
no evidence of CN emission. Based on the sensitiv- 
ity of our observations, we find upper-limit CN and 
H2O production rates of Qcn < 6 x 10 23 mol s _1 
and Qh 2 o < 10 26 mol s _1 , respectively. 

4. Numerical integrations show that P/La Sagra is 
largely dynamically stable, indicating it is likely 
native to its current location in the main belt. The 
surrounding region in orbital element space, how- 
ever, is moderately unstable with two resonances, 
the 13:6 mean-motion resonance with Jupiter and 
the (3,-2,-1) three-body mean-motion resonance 
with Jupiter and Saturn, present nearby. 

5. We find that P/La Sagra is not associated with 
any known asteroid families, though find a small 
cluster of asteroids in its dynamical vicinity. Due 
to the small number of objects in this cluster at the 
current time, however, we cannot confirm whether 
members of this cluster in fact resulted from the 
recent fragmentation of a larger parent asteroid, or 
if they are simply clustered in orbital element space 
by chance. 

Despite our unsuccessful attempt to detect gas emission 
from P /La Sagra, we conclude from the long-lasting and 
steadily evolving morphology of the object's dust cloud 
that its activity is likely to be sublimation-driven, and 
unlikely to be the result of an impulsive impact. As 
such, it is likely to be a genuine comet, making it just 
the fifth MBC discovered (out of a total of six known to 
date). Numerical simulations indicating that P/La Sagra 
is likely to be native to its current location in the main 
belt mean that its physical composition is likely to be rep- 
resentative of other objects in the same region. A search 
for other MBCs in this region of orbital element space 
could prove fruitful. It has been hypothesized, however, 
that present-day MBC activity may require triggering 



Hsieh et al. 


2004 


Hsieh 


2009) 



As 

such, we caution that the lower rate of pot ential a ctivity- 



trigge ring impacts at high inclinations ( Farinella & Davis 
|1992[ ) means that, overall, the rate 01 currently active 
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MBCs in this region is likely to be lower than at lower 
inclinations. 

In any event, the presence of P/La Sagra in a com- 
pletely disparate region of orbital elem ent space fro m the 
other MBCs underscores the caveat bylHsieh (120091) that 



the currently known population is far too small to make 
any conclusions about the abundance and distribution of 
such objects in the asteroid belt. Current and upcoming 
systematic searches of the entire asteroid belt f or MBCs 



now starting to see success in the case of PS1 (Hsieh et 
al.|2011c 2012b ), will be critical for ascertaining the true 



abundance and distribution of icy objects in the inner so 
lar system 
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TABLE 1 
Comet-Like Main-Belt Asteroids 



Name 


Type a 


a b 


e c 


i d 


TjC 


P f 






(A d /A N ) ma J- h 


Discovery Date k 


133P/Elst-Pizarro 


MBC 


3.160 


0.162 


1.386 


3.184 


5.62 


15.49 ±0.05 [1] 


1.9 


0.7 [8] 


1996 Aug 7 [12] 


238P/Read 


MBC 


3.165 


0.253 


1.266 


3.153 


5.63 


19.05 ± 0.05 [2] 


0,1 


21 [2] 


2005 Oct 24 [13] 


176P/LINEAR 


MBC 


3.194 


0.194 


0.238 


3.166 


5.71 


15.10 ±0.05 [1] 


2.0 


0.3 [9] 


2005 Nov 26 [14] 


P/2008 Rl (Garradd) 


MBC 


2.726 


0.342 


15.903 


3.217 


4.50 


20.3 ±0.1 [3] 


0.2 


28 [3] 


2008 Sep 2 [15] 


P/2010 A2 (LINEAR) 


DA 


2.291 


0.124 


5.255 


3.583 


3.47 


21.9 ±0.1 [4] 


0.06 




2010 Jan 6 [16] 


P/2010 R2 (La Sagra) 


MBC 


3.099 


0.154 


21.395 


3.099 


5.46 


17.9 ± 0.2 [5] 


0.7 


60 [5] 


2010 Sep 14 [17] 


(596) Scheila 


DA 


2.928 


0.165 


14.661 


3.209 


5.01 


8.54 [6] 


56.67 


2.2 [10] 


2010 Dec 10 [18] 


P/2006 VW139 


MBC 


3.052 


0.201 


2.438 


3.203 


5.33 


15.9 [7] 


1.7 


1.1 [11] 


2011 Nov 05 [19] 



a Object classification as main-belt comet (MBC) or disrupted asteroid (DA). 

b Osculating semimajor axis, in AU. 

c Osculating eccentricity. 

d Osculating inclination, in degrees. 

e Tisserand parameter. 

f Orbital period, in years. 

s Absolute J?-ba nd magnitude of nucl eus, assuming solar color s, in the IAU H, G system. 

h References: [1] |Hsieh et al.| l |2009a| ; [2] |Hsieh et al.| j2011bl ; [3] jMaclennan fc Hsieh] { 2012\\ [4] |jewitt et al.| |2010[); [5] this work; [6 1 
|Tedesco et~a l. (2004); [71 JPL Small Body Databa se (nttp://ssd. jpl.nasa.gov/sbdb.cgi); |8j Hsieh et al. (2010b); |9| Hs ieh et al^(|2011a|; 

[10] Jcwitt et al. (2011); [111 Hsieh et al. (2012b); [12| h,lst et al. (19y6); |13| Head et ai. (2005); |14| Hsieh et al. (2006); [15] Uarradd et 
[al~1|2UU8| ; [16] |Bi'twhrstle et al.||201U| ; [17] |Nomen et al.| |2U1U| ' |l8j [LarsonppUTU^ [19j [Hsieh et al.||2UTlc| . 



1 Estimated effective radius of nucleus, from same works referenced for Hji, in km. 
j Peak observed dust-to-nucleus scattering surface area ratio. 
k Discovery date of comet-like activity. 
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TABLE 2 
Observation Log 
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2010 Sep 29 


FTN 


N-9 


1.5 


3 


450 


R 


2 


652 


1.790 


13.4 


23.7 


119.9 


238.7 


-11.7 


2010 Oct 05 


Keck 


N-2 


0.9 


2 


240 


BR 


2 


659 


1.824 


14.4 


25.2 


110.7 


238.2 


-11.3 


2010 Oct 05 


Keck 


N-2 


0.9 


4 


3600 


Spec. 


2 


659 


1.824 


14.4 


25.2 


110.7 


238.2 


-11.3 


2010 Oct 19 


UH2.2 


N+ll 


1.1 


3 
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R 


2 


666 


1.922 


16.8 


28.5 


94.5 


237.7 


-9.8 


2010 Nov 26 


UH2.2 


N-9 


1.8 


6 
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R 


2 


698 


2.343 


21.1 


37.8 


71.3 


239.6 


-4.0 


2010 Dec 12 


UH2.2 


N+6 


0.8 


21 


6300 


1! 


2 
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2.535 


21.3 


41.3 


65.7 


241.4 


-1.5 


2011 Dec 31 
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R 


2 
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2.793 


20.5 


45.9 


60.5 


243.9 


1.2 


2011 Feb 03 
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N 


0.8 


2 


360 


BR 


2 


772 


3.210 


17.0 


53.5 


53.4 


248.7 


4.4 


2011 Aug 26 
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N-3 


1.0 


1 
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BR 


3 


066 


3.289 


17.9 


95.1 


270.9 


291.6 


-5.7 


2011 Aug 31 


Gemini 


N+2 


0.8 


6 


1800 


r' 


3 
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3.231 


18.2 


96.0 


270.8 


292.5 


-6.0 


2013 Mar 13 


Aphelion 












3 


570 


2.758 


10.5 


180.0 


76.6 


302.8 


7.5 


2015 Nov 30 


Perihelion 












2 


620 


2.946 


19.3 


0.0 


67.0 


238.0 


-2.8 



a Telescope used (PS1: 1.8-m Pan-STARRS 1 telescope; Dkl.54: Danish 1.54-m telescope; FTN: 2.0-m Faulkes Telescope North; Gemini: 8-m 
Gemini-North telescope; Keck: Keck I 10-m telescope; UH2.2: University of Hawaii 2.2-m telescope; INT: 2.5-m Isaac Newton Telescope) 
b Phase of moon expressed in offset from new moon ("N") in days. 
c Approximate average seeing (FWHM) in arcsec. 
d Number of images. 

e Total effective exposure time in seconds. 
' Heliocentric distance in AU. 
g Geocentric distance in AU. 

h Solar phase angle (Sun - P/La Sagra - Earth) in degrees. 
1 True anomaly in degrees. 

j Position angle of the antisolar vector, as projected in the plane of the sky, in degrees east of north. 

k Position angle of the negative velocity vector, as projected in the plane of the sky, in degrees east of north. 

1 Orbit-plane angle (between observer and object orbit plane as seen from object) in degrees. 



TABLE 3 
Photometry Results 
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Aug 16 
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14.25 ± 0.15 


18.3 ± 


0.1 


14.0 ±0.2 


30 ±10 


40 ± 15 


7±2 


42 ±6 


2010 


Sep 08 


18.26 ± 
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14.20 ±0.15 


18.0 ± 


0.1 


13.9 ±0.2 


30 ± 10 


45 ± 15 


7± 2 


45 ±6 


2010 


Sep 19 (1) 


18.28 ± 
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14.15 ±0.15 


17.9 ± 


0.1 


13.8 ±0.2 


30 ± 10 


45 ± 15 
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46 ±6 


2010 


Sep 19 (2) 


18.23 ± 





02 


14.10 ±0.15 
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0.1 


13.7 ± 0.2 
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45 ± 15 
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48 ±6 


2010 


Sep 29 


18.39 ± 





02 


14.19 ±0.15 


18.0 ± 


0.1 


13.8 ±0.2 


30 ± 10 


45 ± 15 


8±2 


44 ±6 


2010 


Oct 05 


18.34 ± 





03 


14.05 ± 0.20 


18.0 ± 


0.1 


13.7 ± 0.2 


35 ± 10 


50 ±20 


8±3 


49 ± 7 


2010 


Oct 19 


18.47 ± 





03 


13.97 ± 0.20 


18.1 ± 


0.1 


13.6 ±0.2 


35 ± 10 


55 ±20 


9±3 


50 ±8 


2010 


Nov 26 


18.89 ± 





03 


13.80 ± 0.20 


18.4 ± 


0.1 


13.3 ±0.2 


45 ± 15 


65 ±20 


11 ± 4 


49 ±9 


2010 


Dec 12 


18.78 ± 





02 


13.50 ±0.20 


18.4 ± 


0.1 


13.1 ± 0.2 


60 ±20 


85 ±30 


14 ±5 


60 ± 11 


2010 


Dec 31 


18.93 ± 





01 


13.45 ±0.20 


18.2 ± 


0.2 


12.7 ± 0.2 


60 ±20 


85 ±30 


15 ± 5 


57 ± 10 1 


2011 


Feb 03 


19.53 ± 





08 


13.85 ± 0.20 


19.1 ± 


0.1 


13.4 ± 0.2 


40 ± 10 


60 ±20 


10 ±3 


34 ± 6 


2011 


Aug 26 


24.1 ± 





1 


18.1 ±0.3 














0.7 ±0.1 


2011 


Aug 31 


23.4 ± 





1 


17.4 ± 0.3 














1.2 ±0.2 



a Mean apparent magnitude inside a 5''0 photometry aperture. 

b Mean magnitude, normalized to R = A = 1 AU and a = 0°, assuming —0.03 < Gjj < 0.26. 
c Total apparent magnitude inside rectangular aperture enclosing entire comet (nucleus and dust tail) 
d Total magnitude, normalized to R = A = 1 AU and a = 0°, assuming —0.03 < Gr < 0.26. 
e Inferred ratio of scattering surface area of dust to nucleus scattering surface area. 
f Inferred scattering surface area of dust, in 10 6 m 2 , using r^f = 0.6 km. 
s Estimated dust mass, in 10 5 kg, assuming 10 /xm-radius grains and p = 1300 kg m — 3 . 
h Dust contribution (computed using a 5'.'0 photometry aperture), as parameterized by 



A'Hearn et al. 



(1984i, in cm. 



Computed from data in which comet was trailed by 4/1, and as such, represents a lower limit. 
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TABLE 4 
P/La Sagra Cluster Members 



Object 




e b 


sin(i p ) c 




rp e 
Hyap 


18901 (2000 MR 5 ) 


3.0942 


0.1223 


0.3791 


12.1 


683 


106020 (2000 SS294) 
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480 
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0.3824 


15.4 
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P/2010 R2 (La Sagra) 


3.1001 


0.1148 


0.3817 


18.3 


270 



a Proper semimajor axis, in AU 

b Proper eccentricity 

c Sine of proper inclination 

d Absolute magnitude, in V-band 

e Lyapunov time, in kyr 




Semimajor Axis (AU) 

Fig. 1. — Plots of eccentricity (upper panel) and inclination (lower panel) versus semimajor axis showing the distributions in 
orbital element space of main-belt asteroids (black dots), MBCs (red circles), and likely disrupted asteroids (blue circles). Also 
marked with dotted lines are the semimajor axes of Mars (oMars) and Jupiter (aj up ), the semimajor axis of the 2:1 mean-motion 
resonance with Jupiter, and the loci of Mars-crossing orbits and Jupiter-crossing orbits. 
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Fig. 2. — Composite images of P/La Sagra (at the center of each panel) constructed from data obtained on (a) 2010 August 
16 (80 s of effective exposure time on PS1 in r'-band), (b) 2010 September 08 (40 s on PS1 in r'-band), (c) 2010 September 
19 (2400 s on the Danish 1.54 m in .R-band), (d) 2010 September 19 (1500 s on the UH 2.2 m in .R-band), (e) 2010 September 
29 (450 s on FTN in i?-band), (f) 2010 October 05 (240 s on Keck I in i?-band), (g) 2010 October 19 (900 s on the UH 2.2 m 
in ii-band), (h) 2010 November 26 (1800 s on the UH 2.2 m in .R-band), (i) 2010 December 12 (5400 s on the UH 2.2 m in 
.R-band), (j) 2010 December 31 (1200 s on the INT in .R-band), and (k) 2011 February 03 (360 s on Keck I in i?-band). All 
panels are 60" x 60" in size with north (N), east (E), the antisolar direction (—0), and the negative heliocentric velocity vector 
(—v), as projected on the sky, marked. 




Fig. 3. — Composite images of P/La Sagra (point- like object in the lower left of each panel) constructed from data obtained on 
(a) 2011 August 26 (1200 s on Keck I in .R-band) and (b) 2011 August 31 (1800 s on Gemini in r'-band). Panels are 120" x 60", 
with north (N), east (E), the antisolar direction (—0), and the negative heliocentric velocity vector (— v), as projected on the 
sky, marked. 
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Fig. 4. — Orbital position plot of P/La Sagra observations detailed in Table [2] The Sun is shown at the center as a solid dot, with the 
orbits of Mercury, Venus, Earth, Mars, P/La Sagra, and Jupiter (from the center of the plot outwards) shown as black lines. Solid squares 
mark positions where P/La Sagra was observed to be active in 2010 and early 2011, while an open square marks the position where P/La 
Sagra was observed to be inactive in August 2011. Perihelion (P) and aphelion (A) positions are marked with crosses. Observations plotted 
were obtained on (a) 2010 August 16, (b) 2010 September 8-19, (c) 2010 October 5, (d) 2010 November 25-28, (e) 2010 December 12, (f) 
2010 December 31, (g) 2011 February 3, and (h) 2011 August 26-31. 
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Fig. 5. — Comparisons of (a) dust-to-nucleus ratios (by scattering cross-section) a nd (b) total dust s cattering cross-sections, in 10 6 km 2 , 
as func tions of true anomaly m easured for the five known MBCs: 133P (red symb ols; jHsieh et al,|2010b|, 176P (green symbols; |Hsieh et al.| 
|2011a| l, 238P (orange symbols; [Hsieh et al7||2011b[ | , P/Garradd (purple symbols; [Jewitt et al.|ijOUy| lVlaclennan & Hsich 2012|, and P/La 
sagra (blue symbols; this workj! Ubservations wnere activity was detected either visually or from photometry are marked with circular 
symbols, while observations where no activity was detected are marked with squares. For points where no error bars are visible, the amount 
of uncertainty is equal to or less than the size of the plotted symbol. 
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Fig. 6. — Plots of the sky-background-subtracted and flux-calibrated spectrum of P/La Sagra (solid black line), and the scaled spectrum 
of the solar analog, the G2V star HD 91163 (blue dashed line). The shaded region indicates the wavelength region where the CN emission 
band is expected. We note that the spectrum of P/La Sagra closely resembles that of HD 91163 and thus no emission features are observed 
within the uncertainties of our measurements. 




Fig. 7. — The spectrum of P/La Sagra with the underlying solar continuum removed. The three red error bars show the 1-cr uncertainties 
in the three wavelength regions (3760-3830A, 3830-3900A, and 3900-3970A). 
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Fig. 8. — Plots of semimajor axis versus eccentricity (left) and inclination (right) showing initial osculating elements of test p articles in 
three sets of 100 Gaussian-distributed particles each of which are subjected to a 100 Myr dynamical integration (Section |3.4.1[ l. Each set 
of test particles is centered on the current osculating orbital elements of P/La Sagra (red circle) and characterized by a a value equal to 
the object's JPL-tabulated uncertainties. Particles ejected in less than 20 Myr, between 20 Myr and 50 Myr, and between 50 Myr and 100 
Myr are plotted with orange X symbols, yellow asterisk symbols, and green triangles, respectively, while particles that are not ejected after 
the 100 Myr test period used in our simulations are marked with blue squares. 
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Fig. 9. — Same as Figure [8] for three sets of 100 Gaussian-distributed particles each, where each set is characterized by a a value equal 
to 100X P/La Sagra's JPL-tabulated uncertainties. The locations of the (3,-2,-1) three-body mean-motion resonance (with Jupiter and 
Saturn) at a = 3.0790 AU and the 13:6 mean-motion resonance (with Jupiter) at a = 3.1063 AU are plotted with dashed lines. 
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Fig. 10. — Plot of P/La Sagra's semimaj or axis (AU) as a function of time (Myr), marked by a black solid line, as simulated in the 
numerical integrations described in Section |3.4.1| where the semimajor axes of the 13:6 mean-motion resonance with Jupiter and the 
(3,-2,-1) three-body mean-motion resonance with Jupiter and Saturn arc marked with red dashed lines. 
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Fig. 11. — Synthetic orbital element plots of eccentricity (upper panel) and sine of inclination (lower p anel) versus semim ajor axis of 
outer main-belt asteroids (small black dots), where the central bodies of families and clumps identified by Gil-Hutton (2006) a re ma rked 
with blue and green triangles, respectively, and central bodies of families, clumps, and clusters identified by JNovakovic et alH |20TTJ are 
marked with blue, green, and orange circles, respectively. The location of P/La Sagra in synthetic proper element space is plotted with a 
red circle. 
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Fig. 12. — Number of objects associated with the P/La Sagra cluster. The six members determined to belong to this cluster are all 
linked within 95 m s^ 1 , with no more linkages found to other asteroids until 125 m s^ 1 when the group merges with the local background 
population. Our choice for the nominal cutoff distance for the cluster of 110 m s — 1 is marked. 
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Fig. 13. — Average differences (AQ) of the longitudes of the ascending nodes for members of the P/La Sagra cluster. Three significant 
clusterings at 4.9, 6.2 and 18.2 Myr in the past are marked with arrows. 



